Astronomy & Astrophysics manuscript no. aa 


©ESO 2012 


October 16, 2012 





Full-disk nonlinear force-free field extrapolation of SDO/HMI and 

SOLIS/VSM magnetograms 

Tilaye Tadesse^'^, T. Wiegelmann^, B. Inhester^, P. MacNeice'*, A. Pevtsov^, and X. Sun^ 

' Department of Physics, Drexel University, Philadelpliia, PA 19104-2875, USA 

e-mail: tasf awOeinstein . physics . drexel . edu 
' Addis Ababa University, Institute of Geophysics, Space Science, and Astronomy, Po.Box 1 176, Addis Ababa, Ethiopia 

e-mail: tilaye . tadesse@gmail . com 
^ Max-Planck-Institut fiir Sonnensystemforschung, Max-Planck-Strasse 2, 37191 Katlenburg-Lindau, Germany 

NASA/ GSFC, Greenbelt, MD, U.S.A. 
^ National Solar Observatory, Sunspot, NM 88349, U.S.A. 

* W.W. Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA 
Preprint online version: October 16, 2012 

ABSTRACT 

Context. The magnetic field configuration is essential to understand solar explosive phenomena such as flares and coronal mass 
ejections. To overcome the unavailability of coronal magnetic field measurements, photospheric magnetic field vector data can be 
used to reconstruct the coronal field. A complication of this approach is that the measured photospheric magnetic field is not force- 
free and that one has to apply a preprocessing routine in order to achieve boundary conditions suitable for the force-free modelling. 
Furthermore the nonlinear force-free extrapolation code should take uncertainties into account in the photospheric field data which 
occur due to noise, incomplete inversions or azimuth ambiguity removing techniques. 

Aims. Extrapolation codes in cartesian geometry for modelling the magnetic field in the corona do not take the curvature of the 
Sun's surface into account and can only be applied to relatively small areas, e.g., a single active region. Here we apply a method for 
nonlinear force-free coronal magnetic field modelling and preprocessing of photospheric vector magnetograms in spherical geometry 
using the optimization procedure to full disk vector magnetograms. We compare the analysis of the photospheric magnetic field 
and subsequent force-free modeling based on full-disk vector maps from Helioseismic and Magnetic Imager (HMl) on board solar 
dynamics observatory (SDO) and Vector Spectromagnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun 
(SOLIS). 

Methods. We use Helioseismic and Magnetic Imager and Vector Spectromagnetograph photospheric magnetic field measurements to 
model the force-free coronal field above multiple solar active regions, assuming magnetic forces to dominate. We solve the nonlinear 
force-free field equations by minimizing a functional in spherical coordinates over a full disk excluding the poles. After searching 
for the optimum modeling parameters for the particular data sets, we compare the resulting nonlinear force-free model fields. We 
compare quantities like the total magnetic energy content and free magnetic energy , the longitudinal distribution of the magnetic 
pressure and surface electric current density using our spherical geometry extrapolation code. 

Results. The magnetic field lines obtained from nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager and 
Vector Spectromagnetograph data have good agreement. However, the nonlinear force-free extrapolation based on Helioseismic and 
Magnetic Imager data have more contents of total magnetic energy, free magnetic energy, the longitudinal distribution of the magnetic 
pressure and surface electric current density compared to the one from Vector Spectromagnetograph data. 
Conclusions. 
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1. Introduction 

The magnetic fields configuration is essential for us to under- 
stand solar explosive phenomena such as flares and coronal mass 
ejections. The corona has been the subject of extensive modeling 
for decades, but these efforts have been hampered by our lim- 
ited abihty to determine the corona's three-dimensional struc- 



ture (Schrijver & Title 2011 [Sandman & Aschwanden 2011 



Since the corona is optically thin, direct measurements of these 
magnetic fields are very difficult to implement, and the present 
observations for the magnetic fields based on the spectropolaii- 
metric method (the Zeeman and the Hanle effects) are limited to 
low layers of solar atmosphere (photosphere and chromosphere). 
The problem of measuring the coronal field and its embedded 
electrical currents thus leads us to use numerical modelling to in- 



fer the field strength in the higher layers of the solar atmosphere 
from the measured photospheric field. Due to the low value of 
the plasma /3 (the ratio of gas pressure to magnetic pressure), 
the solar corona is magnetically dominated ( Gary|200T] l. To de- 
scribe the equilibrium structure of the static coronal magnetic 
field when non-magnetic forces are negligible, the force-free as- 
sumption is appropriate: 

(VxB)xB = (1) 

V ■ B = (2) 

subject to the boundary condition 

B = Bobs on photosphere (3) 

where B is the magnetic field and Bobs is measured vector field 
on the photosphere. Equation ([T]i states that the Lorentz force 
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vanishes (as a consequence of J || B, where J is the elec- 
tric current density) and Equation (|2| describes the absence of 
magnetic monopoles. Based on the above assumption, the coro- 
nal magnetic field is modelled with nonlinear force-free field 
(NLFFF) extrapolation (Inhester & Wiegelmann'12006', 'Valoril 



et al. 



2005 1 Wiegebnann 



2004 Wheatland 2004; Wheatland & 



Regnier||2009 ', Tadess e et al.||2009| [Wheatland & LekapOll 



Amari & Aly 2010, Wiegeknann et al. 2012 Jiang & Feng 
2012| l. From a mathematical point of view appropriate boundary 



condition for force-free modehng are the vertical magnetic field 
B„ and the vertical current J„ prescribed only for one polarity of 
B„ (Ama ri etal. 1 19 97 1999 /2006). A direct use of t hese bound- 



ary conditions is implemented in Grad-Rubin codes ( Amari et al 
1999]! [Wheatland & Regnier] ( p009) and [Wheatland & Leka 



( 201 1 1 implemented the use of Z?,^ and Z?„ solution together with 



an error approximation to derive consistent solutions. Using the 
three components of B as boundary condition requires consistent 



magnetograms, as outlined in Aly ( 1989) . We use preprocessing 
and relaxation of the boundary condition to derive these consis- 
tent data on the boundary. 

As an alternative to real measurement, nonlinear force-free 
field (NLFFF) models are thought to be viable tools for investi- 
gating the structure, dynamics, and evolution of the coronae of 
solar active regions. It has been found that NLFFF models are 
successful in application to analytic test cases (Schrijv er et al.[ 
2006 Metcalf et al. 2008 1, but they are less successful in ap- 
phcation to real solar data. However, NLFFF models have been 
adopted to study various magnetic field structures and proper - 
ties in the solar atmosphere. For instanc e, [Regnier et al.[( 2002i; 

('20D4^;'Canou et al. (2009); Canou & Amari 
(2010) ; Valori et al. (2012) have substantially 



[Regnier & Amari 



(2010) ; ,Guo etal. 



studied various magnetic field structures and properties using 
their respective NLFFF model codes. Different NLFFF models 
have been found to have markedly different field line configu- 
rations and to provide widely varying estimates of the magnetic 
free energy in the coronal volume, when applied to solar data 
( [DeRosa et al.[[2009| ). The main reasons for that problem are 
(1) the forces acting on the field within the photosphere, (2) the 
uncertainties on vector-field measurements, particularly on the 
transverse component, and (3) the large domain that needs to be 
modelled to capture the connections of an active region to its 
surroundings) Tade sse et al.|201 1 2012a| l. In this study, we have 
considered those three points explicitly into account. However, 
caution must still be needed while assessing results from this 
modeling. This is because many aspects of the specific approach 
to modeling used in this work, such as the use of preprocessed 
boundary data, the missing boundary data, and the departure of 
the model fields from the observed boundary fields may influ- 
ence the results. 

In this work, we use full-disk SDO/HMI and SOLIS/VSM 
photospheric magnetic field measurements to model the NLFFF 
coronal field above multiple solar active regions. Comparison of 
vector magnetograms for one particular active region observed 
with two different instruments from SOLIS and HMI and their 
corresponding force-free models have been studied by [Thalmann[ 
|et al.[ ( [20T2] submitted to AJ) in Cartesian coordinates. We use a 
larger computational domain which accommodates most of the 
connectivity within the coronal region. We use a spherical ver- 
sion of the optimization procedure that has been implemented 
in Tadesse et al. ( 201 l| l. We compare quantities like the total 
magnetic energy content and free magnetic energy and the lon- 
gitudinal distribution of the magnetic pressure in the HMI and 
VSM-based model volumes in spherical geometry. We relate the 
appearing differences to the photospheric quantities such as the 



magnetic fluxes and electric currents but also show the extent of 
agreement of NLFFF extrapolations from different data sources. 



2. Instrumentation and data set 

2A. Solar Dynamics Observatory(SDO) 
Magnetic Imager(HMI) 



Helioseismic and 



The Helioseismic and Magnetic Imager (HMI; Schou et al. 
[201 2p is part of the Solar Dynamics Observatory (SDO) and ob- 
serves the full Sun at six wavelengths and full Stokes profile in 
the Fe I 617.3 nm spectral line. HMI consists of a refracting 
telescope, a polarization selector, an image stabilization system, 
a narrow band tunable filter and two 4096 pixel CCD cameras 
with mechanical shutters and control electronics. Photospheric 
line-of-sight LOS and vector magnetograms are retrieved from 
filtergrams with a plate scale of 0.5 arc-second. From filter- 
grams averaged over about ten minutes. Stokes parameters are 
derived and inverted using the Milne-Eddington (ME) inversion 
algorithm of Borrero et al. (2011) (the filling factor is held at 
unity). Within automatically identified regions of strong mag- 
netic fluxes (Turmon et al. [2010) 1, the full disk inversion data 
are from the second HMI vector data release (JSOC data series 
hmi.ME_720s_el5wl332). The 180-degree azimuthal ambiguity 
in the strong field region is resolved using the Minimum Energy 
Algorithm (M etcalf 199 4; Metcalf et al. 2006; Leka et al. 2009) , 
taken from the AR patches in the second release also (data series 
hmi.B_720s_el5wl332_cutout). For the weak field region where 
noise dominates, we adopt a radial-acute angle method to re- 
solve the azimuthal ambiguity. The weak field region is defined 
as where field strength is below 200 G at disk center, 400 G on 
the limb, and varies linearly in between. The noise level is » lOG 
and a: lOOG for the longitudinal and transverse magnetic field, 
respectively. 

2.2. Synoptic Optical Long-term Investigations of the 
Sun(SOLIS) - Vector-SpectroMagnetograph(VSM 



The Vector Spectromagnetograph (VSM; see [Jones et al.|2002| l 
is part of the Synoptic Optical Long-term Investigations of the 
Sun (SOLIS) synoptic facility (SOLIS; see |Keller et al.||2003] l. 
VSM is a full disk Stokes polarimeter As part of daily synoptic 
observations, it takes four different observations in three spec- 
tral lines: Stokes /(intensity), V (circular polarization), Q, and 
U (linear polarization) in photospheric spectral lines Fe 1 630. 15 
nm and Fe I 630.25 nm , Stokes / and V in Fe I 630.15 nm and 
Fe I 630.25 nm, similar observations in chromospheric spectral 
line Ca II 854.2 nm, and Stokes / in the He 1 1083.0 nm line and 
the near by Si I spectral line. Observations of /, Q, U, and V are 
used to construct full disk vector magnetograms, while I-V ob- 
servations are employed to create separate full disk longitudinal 
magnetograms in the photosphere and the chromosphere. The 
vector data are provided with a plate scale of one arc-second. 
The lower limits for the noise levels are a few Gauss in the lon- 
gitudinal and 70G in the transverse field measurements. 

Quick-look (QL) vector magnetograms were created based 
on an algorithm by Auer et al. ( 1977 1. Beginning January 2012, 
QL vector magnetograms are created using weak-field approx- 
imation ^Ronan et al. |1987[ l. The algorithm uses the Milne- 
Eddington model of solar atmosphere, which assumes that the 
magnetic field is uniform (no gradients) through the layer of 
spectral line formation (Unno 1956). It also assumes symmetric 
line profiles, disregards magneto-optical effects (e.g., Faraday 
rotation), and does not distinguish the contributions of mag- 
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netic and non-magnetic components in spectral line profiles (i.e., 
magnetic filling factor is set to unity). A complete inversion of 
the spectral data is performed later using a technique developed 
by|Skumanic h & Lites| ( fl987l ). This latter inversion (called ME 
magnetogram) also employs Milne-Eddington model of atmo- 
sphere, but solves for magneto-optical effects and determines the 
magnetic filling factor i.e., (the fractional contribution of mag- 
netic and non-magnetic components to each pixel). The ME in- 
version is only performed for pixels with spectral line profiles 
above the noise level. For pixels below the polarimetric noise 
threshold, the magnetic field parameters are set to zero. 

From the measurements, the azimuths of transverse mag- 
netic field can be determined with 180-degree ambiguity. This 
ambiguity is resolved using the non-potential field calculation 
(NPFC; see |Georgoulis|2005| ). The NPFC method was selected 
on the basis of a comparative investigation of several methods 
for 180-degree ambiguity resolution (Met calf et al.| 2006). Both 
QL and ME magnetograms can be used for potential and/or 
force-free field extrapolation. However, in strong fields inside 
sunspots, the QL field strengths may exhibit an erroneous de- 
crease inside the sunspot umbra due to so-called magnetic satu- 
ration. For this study, we choose to use fully inverted ME mag- 
netograms. 

3. Method 

Photospheric field measurements are often subject to measure- 
ment errors. In addition to this, there are finite non-magnetic 
forces which make the data inconsistent as a boundary for a 
force-free field in the corona. In order to deal with these un- 
certainties, one has to: 1.) preprocess the surface measurements 
in order to make them compatible with a force-free field and 2.) 
keep a balance between the force-free constraint and deviation 
from the photospheric field measurements. Both methods con- 
tain free parameters, which have to be optimized for use with 
data from SOLIS/VSM and SDO/HMI. 

3. 7. Preprocessing of HMI and VSM data 

To serve as suitable lower boundary condition for a force-free 
modeling, vector magnetograms have to be approximately flux 
balanced and on average a net tangential force acting on the 
boundary and shear stresses along axes lying on the boundary 
have to reduce to zero. We use dimensionless parameter s, efiux, 
eforce and etorque, to quanti fy such properties(Wiegelmann et al. 
[200^ [Tadesse et al.|2009||Aly||1989tpolodensky,1969| . Even 
if we choose a sufficiently flux balanced region (ejiux), we find 



that the force-free conditions e 



-force 



<K 1 and e, 



lorque 



<K 1 are not 



usually fulfilled for measured vector magnetograms. In order to 
fulfill those conditions, we use preprocessing method as imple- 
mented in Wiegelmann et al.| (2006i. The preprocessing scheme 
of |Tadesse et al. (2009) involves minimizing a two-dimensional 
functional of quadratic form in spherical geometry similar to 

B — argmin(Lp), 

Lp = yUiLi -H 112L2 + n-iL^ + haU, (4) 

where B is the preprocessed surface magnetic field from the in- 
put observed field Bots- Each of the constraints L„ is weighted by 
an as yet undetermined factor /i„. The first term (n = 1) corre- 
sponds to the force-balance condition, the next (n = 2) to the 
torque-free condition, and the last term (n - 4) controls the 
smoothing. The explicit form of Li, L2, and L4 can be found 
in |Tadesse et a l. (2009). The term (n - 3) controls the difference 
between measured and preprocessed vector fields.. 



3.2. Optimization principle 

We solve the force-free equations ([T) and ([2| by optimization 
principle, as propose d by [Wheatland et aT] ( |2000| l and general- 
ized by Wiegelmann|(|2004|l for cartesian geometry. The method 
minimizes a joint measure of the normalized Lorentz forces and 
the divergence of the field throughout the volume of interest, 
V. Throughout this minimization, the photospheric boundary of 
the model field B is matched exactly to the observed Bohs and 
possibly preprocessed magnetogram values B. Here, we use the 
optimization approach for functional (L^) in spherical geome- 
try ( Wiegelmann|2007 Tadesse et al.|2009p along with the new 
method, which instead of an exact match enforces a minimal de- 
viation between the photospheric boundary of the model field B 
and the magnetogram field Bobs by adding an appropriate sur- 
face integ ral term Lp/,0,0 ( Wiegelmann & Inhester|2010 Tadesse 
et al.|201 1| ). These terms are given by 



B — argmm(Li^) 



La, - Lf + + vLphoig 



(5) 



-'photo 



Jv 



"^|(V xB)x B^P- smedrd0d4> 



■Bfr^ sm0drd0d(f> 



Bobs) ■ W(0, <p)-{B- Bobs)r^ sin 0d0d(p 



where Lf and measure how well the force-free Eqs. (fill 
and divergence-free (j2]i conditions are fulfilled, respectively, 
and both <wy (r, 0, <p) and ciJdir, 0, (p) are weighting functions. The 
weighting functions coj and ajj in Lf and in Eq. (|5]l are cho- 
sen to be unity within the inner physical domain V and dechne 
with a cosine profile in the buffer boundary region ( [Wiegelmann[ 
|2004; ,Tadesse et al.|[2009 2012b). They reach a zero value at 
the boundary of the outer volume V. The distance between the 
boundaries of V and V is chosen to be nc/ = 10 grid points wide. 
The third integral, Lphoro, is the surface integral over the photo- 
sphere which allows us to relax the field on the photosphere to- 
wards force-free solution without too much deviation from the 
original surface field data. 

W(0, 0) is a space-dependent diagonal matrix the element of 
which are inverse proportional to the estimated squared mea- 
surement error of the respective field component. In principle 
one could compute W from the measurement noise and er- 
rors obtained from the inversion of measured Stokes profiles 
to field components. Until these quantities become available, 
a reasonable assumption is that the magnetic field is measured 
in strong field regions more accurately than in the weak field 
and that the error in the photospheric transverse field is at 
least one order of magnitude higher as the line-of-sight com- 
ponent. Appropriate choices to optimize v and W for use with 
SDO/HMI(Wiegelmann et al.|[2012|) and SOLIS/VSM( Tadesse 



et al. 201 l| l magnetograms have been investigated. For a de- 



tailed description of the current code implementation, we refer 
to Wiegelmann & Inhester ( 2010 1 and Tadesse et al. ( 201 1 1. 



4. Results 

Within this work, we use the full disk data from SOLIS/VSM 
and SDO/HMI instruments obsrved on November 09 2011 
around 17:45UT. During this observation there were four ac- 
tive regions (ARs 11338, 11339, 11341 and 11342) along with 
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Fig. 1. Magnetic vector maps of VSM and HMI on part of the lower boundary. The color coding shows Br on the photosphere and 
the white arrow indicates the transverse components of the field. The vertical and horizontal axes show latitude, Q and longitude, 
on the photosphere respectively. 



other smaller sunspots spreading on the disk. To accommodate 
the connectivity between those ARs and their surroundings, we 
adopt a non uniform spherical grid r, 6, with rir - 225, 
rig - 375, = 425 grid points in the direction of radius, 
latitude, and longitude, respectively, with the field of view of 
[r„i„ = IRq : r^ax = 2/?0]x[6'„,„ = -50° : fl^ax = 5O°]x[0„in = 
90° : 0max = 270°]. Given the twice as large plate scale of VSM, 
we bin the HMI vector maps to the resolution of VSM in or- 
der to compare the photospheric magnetic field and subsequent 
force free modeling. To deal with vector magnetogram data be- 
ing inconsistent with the force-free assumption, we use a pre- 
processing routine in spherical geometry, which derives suitable 
boundary conditions for force-free modeling from the measured 
photospheric data. Applying this procedure to both SDO/HMI 
and SOLIS/VSM reduces e/orce and e,orq„e further significantly. 
The two quantities are very well below unity after preprocess- 
ing, which gives us some confidence that the data might serve 
as suitable boundary condition for a force-free modeling. Doing 
this, we do not intend to suppress the existing forces in the pho- 
tosphere. Instead, we try to approximate the magnetic field at 
a chromospheric level where magnetic forces are expected to 
be much smaller than in the layers below. Both vector magne- 
tograms are almost flux balanced and the field of view was large 
enough to cover the full-disk. The unsigned magnetic flux of lon- 
gitudinal surface magnetic field from HMI is 1 .57 times that of 



VSM magnetogram. This is in agreement with recent compara- 
tive study by Pietarila et al. (2012 ), who found that the factor to 



convert SOLIS/VSM to SDO/HMI increases with flux density 
from about 1 (weak fields) to about 1.5 (strong fields) for the 
line-of-sight full disk magnetograms. HMI inverts weak field re- 
gions, however, for VSM zeros are assigned to pixels where the 
measured polarization signal is too weak to perform a reliable 
inversion. Disregarding these zero-pixels, about 20% of the total 
number of pixels in the HMI and VSM full disk vector maps are 
remaining for comparison. HMI is found to detect most trans- 
verse field. 

We used a standard preprocessing parameter set ^\ - ^2 - ^ 
and //3 = 0.001, which are similar to the values calculated from 
vector data used in previous studies( Wiegelmann et al. 2012 1 for 
HMI data in Cartesian coordinates. Table [T] lists the values of 
dimensionless parameters for the used HMI and VSM data-sets. 
In this study, we have found that the optimal value for smoothing 
parameter is 114 - 0.05 for full-disk HMI data. These parameters 
control the amount of force-freeness, torque-freeness, nearness 
to the actually observed data and smoothing, respectively. As 



the result of parameter study, Tadesse et al. (2011 1 have found 
A*! - ^J^2 - L Afs = 0.03 and /14 - 0.45 as optimal for full-disk 
VSM data. 

The preprocessing influences the structure of the magnetic 
vector data. It does not only smooths B, (transverse field) but 
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Table 1. Flux-balance, force and torque free parameters of 
SOLIS/VSM and SDO/MHI full disk magnetograms. 



Data set 






^torque 


HMI observed 


-0.0621 


0.1305 


0.1773 


HMI preprocessed 


-0.0313 


0.0001 


0.0002 


SOLIS observed 


-0.0857 


0.4571 


0.2947 


SOLIS preprocessed 


-0.0460 


0.0015 


0.0007 



also alters its values in order to reduce the net force and torque. 
The change in B, is more pronounced than the radial compo- 
nent Br (radial field) since B, is measured with lower accuracy 
than the longitudinal magnetic field. Figure [T] shows the prepro- 
cessed and observed surface vector magnetic field obtained from 
SDO/HMI and SOLISA'SM magnetograms. To identify the sim- 
ilarity of vector components from HMI and VSM on the bottom 
surface, we calculate their pixel-wise correlations before and af- 
ter preprocessing. The correlation were calculated from 



C 



1/2' 



(6) 



where v-, and m,- are the vectors at each grid point / on the bottom 
surface. If the vector fields are identical, then C„pc = 1; if fi -L 
, then C„ec = 0. Table |2] shows the correlation {C„c) of the 2D 
surface magnetic field vectors of observed and preprocessed data 
from HMI and VSM for the radial and transverse components. 
The vector correlation between B, in the preprocessed HMI and 
VSM surface vector maps is clearly more closer to unity than the 
corresponding surface vector maps without preprocessing. There 
is no such difference in correlations between By before and af- 
ter preprocessing. This is to be expected since the preprocessing 
scheme only smooths the longitudinal field while it smooths and 
alters the transverse field. The mean value of the changes due 
to preprocessing in the longitudinal field is lO^^^G and for the 
transverse field on the order of 10 G, i. e., well within the mea- 
surement uncertainty of the HMI and VSM. 



Before we perform nonlinear force-free extrapolations we use 
the preprocessed radial component B, of the VSM and HMI- 
data to compute the corresponding potential fields using spher- 
ical harmonic expansion for initializing our code. We imple- 
ment the new term Lphoto in Eq. dSll to work with boundary 
data of different noise levels and qualities or even neglect some 
data points completely. SOLIS/VSM provides full-disk vector- 
magnetograms, but for some individual pixels the inversion from 
line profiles to field values may not have been successfully in- 
verted and field data there will be missing for these pixels. Since 
the old code without the Lphoto term requires complete boundary 
information, it cannot be applied to this set of SOLIS/VSM data. 
In our new code, these data gaps are treated by setting W — ^ for 
these pixels in Eq. (j5]l ( [Wiegelmann & InhesterpOlO Tadesse 
|et al.|[20TT] l. For those pixels, for which Bobs was successfully 



inverted, we allow deviations between the model field B and the 
input fields observed Bobs surface field using Eq. (j5]l, so that the 
model field can be iterated closer to a force-free solution even if 
the observations are inconsistent. 

The improved optimization scheme allows us to relax the 
magnetic field also on the lower boundary. The relaxation of the 
lower boundary introduces a further modification of the vector 
data, in addition to that by the preprocessing applied before. The 
mean modification of the longitudinal field due to the relaxation 




la) SOLIS/VSM 




Fig. 2. Mask function for magnetic vector field distribution on 
full disk from (a) VSM and (b) HMI. The vertical and horizon- 
tal axes show latitude, Q and longitude, on the photosphere 
respectively. 



of the lower boundary is 10 "^G and absolute values are on the 
order of 1 G. The mean changes of the transverse field are on the 
order of lOG and absolute values can be several 100 G. Given the 
noise levels of HMI and VSM measurements of the longitudinal 

lOG and a few G, respectively) and transverse field (a; lOOG 
and > 70G, respectively), the modifications are on the order of 
the measurement error. 

For nonlinear force-free fields we minimize the functional 
Eq. (|5]). In order to control the speed with which the lower 
boundary is injected during the extrapolation, we vary the 
Langrangian multiplier v. Unless an exact error computation be- 
comes available from inversion and ambiguity removal of the 
photospheric magnetic field vector, a reasonable assumption is 
that the field is measured more accurately in strong field regions 
and one can carry out computations with the mask oc B, and 
oc B]. We choose a mask function of W = {Btlmax{B,)f , which 
gives more weight to strong field regions than to weak ones as 
investigated in Wiegelmann et al. ( 2012[ ). Figure |2] shows sur- 
face distribution of mask function W for VSM and HMI fall-disk 
data. For strong field regions W is close to unity and decline to 
zero in weaker field regions. We vary the Langrangian multiplier 
V between 0.1 and 0.0001 to investigate the optimal parameter 
for HMI full-disk data. To evaluate how well the force-free and 
divergence-free condition are satisfied for different Langrangian 
multiplier v, we monitor a number of expressions, such as Lj, Ld 
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Table 2. The correlations between the components of surface 
fields from HMI and VSM data. 





V 


u 


Cvec 


No preprocessing 


{BHMl)r 


(BvSM)r 


0.947 


No preprocessing 


(Bhmi)i 


(Bvsm)i 


0.893 


Preprocessed 


(BHMl)r 


{Bvsm),- 


0.965 


Preprocessed 


(Bhmi)i 


(Bvsm)i 


0.951 



Table 3. Evaluation of force-free field models from preprocessed 
HMI data. The first column names the model and in column 2 
shows the used Langrangian multipliers. Columns 3-5 show dif- 
ferent force-free consistency evaluations. Column 6 shows the 
ratio of NLFFF energy density to the corresponding potential 
energy density and column 7 the computing time. 



V 


^/ 


Ld 


sin '(o"i) 


^/ ^pot 


Time 


0.1 


21.7 


13.4 


25.8° 


1.06 


2h:17min 


0.05 


19.8 


10.7 


18.1° 


1.12 


3h:31 min 


0.001 


2.9 


1.5 


4.8° 


1.22 


4h:39min 


0.005 


5.2 


3.9 


8.9° 


1.23 


llh:47min 


0.0001 


7.7 


4.3 


10.2° 


1.26 


48h:53min 



Table 4. The magnetic energy associated with extrapolated 
NLFFF field configurations from full disk SDO/HMI and 
SOLIS/VSM data. 



Model 


E„iff{W"erg) 


£ftee(10^^erg) 


SOLIS/VSM 


8.609 


1.375 


SDO/HMI 


8.913 


1.607 



and 



Bi 



(7) 



where crj is the sine of the current weighted average angle be- 
tween the magnetic field B and electric current /. 

For a sufficient small Lagrangian multiplier v - 0.001 we 
found that the resulting coronal fields are force and divergence 
free compared to other values as shown in Table[3] The weighted 
angle between the magnetic field and electric current is about 5° 
for V = 0.001 . Injecting the boundary faster by choosing a higher 
Lagrangian multiplier (v - 0.1) speeds up the computation, but 
the residual forces are higher and current and field are not well 
aligned as investigated by Wiegebnann et al. (2012 1 for single 
AR. 

To understand the physics of solar flares, including the lo- 
cal reorganization of the magnetic field and the acceleration of 
energetic particles, one has to estimate the free magnetic en- 



ergy available for these phenomena (Regnier & Priest 2007a 
Aschwanden|2008 Schrijver|2009 1. This is the free energy that 



can be converted into kinetic and thermal energy. From the en- 
ergy budget and the obs erved magnet ic activity in the active re- 
gion, Regnier & Priest ( 2007b) and |Thalmann et al. ( 2008) in- 
vestigated the free energy above the minimum-energy state for 
the flare process. We estimate the free magnetic energy to be 
the difference between the extrapolated force-free fields and the 
potential field with the same normal boundary conditions in the 
photosphere. We therefore estimate the upper limit to the free 
magnetic energy associated with coronal currents of the form 



-B 



por)'"^ sin0drd6d(p. 



(8) 



where Bpot and B,,/// represent the potential and NLFFF mag- 
netic field,respectively. The magnetic energy densities asso- 
ciated with the potential field configurations from SDO/HMI 
and SOLIS/VSM data are found to be 7.306 x \(fi^erg and 
7.234 X IQ^^erg, respectively. This has to be expected as the un- 
signed magnetic flux of longitudinal surface magnetic field from 
HMI is greater than that of VSM magnetogram. The magnetic 
energy of NLFFF obtained from HMI data is greater that the 
one obtained from VSM data as shown in Table |4] This is due 
to the fact that HMI data has more longitudinal unsigned mag- 
netic flux and detects more transverse field than VSM. To study 
the influence of the use of preprocessed boundary data from the 
observed boundary fields on the estimation of free-magnetic en- 
ergy , we have computed the magnetic energy associated with 
the potential field and NLFFF configurations from the origi- 
nal SDO/HMI data without preprocessing and with preprocess- 
ing. The case for SOLIS/VSM has been studied by Tadesse] 
et al. ( 2012a| l. As preprocessing procedure filters out small scale 



surface field fluctuations, the magnetic energy associated with 
NLFFF obtained from preprocessed SDO/HMI boundary data is 
smaller than the one without preprocessing. Obviously, the po- 
tential field energies of boundary data with and without prepro- 
cessing are close in value, since the potential field calculation 
makes use of the radial magnetic field component which is not 
affected too much by preprocessing procedure. The computed 
magnetic energy from SDO/HMI original data without prepro- 
cessing is about 9.067 x lO-'^erg, which is about 1.7% higher 
than the one obtained from preprocessed and modified observa- 
tional HMI boundary data. However, this energy does not cor- 
respond to the nonlinear force-free magnetic field solution since 
the original boundary data without preprocessing is not a con- 
sistent boundary condition for NLFFF modeling! Tadesse et al. 
|2012a| l. 

We investigate the magnetic field configurations of the VSM 
and HMI models by comparing the vector components. We cal- 
culate the vector correlation (using Equation (j6]l in the com- 
putational volume) of the potential fields and the NLFFF fields. 
The average vector correlation between the potential fields based 
on the HMI and VSM data is 0.97. The average vector correla- 
tion between the NLFFF fields of HMI and VSM data is 0.94. 
Figure [3] a. and b. show the surface radial magnetic field com- 
ponent observed by HMI instrument on November 09 2011 and 
the corresponding AIA (Atmospheric Imaging Assembly) image 
in 171 A, respectively. The magnetic field lines obtained from 
nonlinear force-free extrapolation based on HMI and VSM data 
have good correlations as shown in Figure |3] c. and d., with 
the foot points of the field lines from the two magnetograms 
are identical. However, there are some differences. For example, 
extrapolated field lines from SDO/HMI magnetogram (Figure 
3d) do not show transequatorial loops connecting trailing polar- 
ity of NOAA AR 11339 (west of central meridian in northern 
hemisphere) and trailing polarity of AR11338 (southern hemi- 
sphere). This transequatorial loop is well represented by NLFFF 
extrapolation based on SOLIS/VSM. This difference can be at- 
tributed to presence of a patch of weak fields between two ac- 
tive regions (in SDO/HMI data). With this weak field patch, 
SDO/HMI model tends to close field lines originating in trail- 
ing polarity of AR11339, while SOLIS/VSM model extends 
them to ARl 1338. Both extrapolations indicate loops connecting 
AR 11339 and ARl 1342 (East of central meridian in Northern 
hemisphere). Although SDO/AIA image (Figure 3b) does not 
show coronal loops connecting these two active region, such 
loops are clearly visible in images taken by X-ray Telescope 
on Hinode. These loops appear to fit better field lines from 
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(c) NLFFF solulion from SOLIS/V SM (d ) NLFFF solution from SDO/HMI 



Fig. 3. a) SDO/HMI and b) AIA images and their respective selected magnetic field lines plots reconstructed from c) SOLIS and 
d) HMI magnetograms using nonlinear force-free modelling. The color coding shows on the photosphere. Yellow field lines 
represent closed field lines, while field lines changing in color from yellow to brown (from bottom to the top) represent the open 
ones. The gray area indicates the region where magnetic field values are close to zero. 



SOLIS/VSM model. Despite a relatively good visual agreement 
in extrapolated fields, the models show some notable disagree- 
ment in derived magnetic energy. Thus, for example, the esti- 
mated free magnetic energy obtained from SDO/HMI is 14.4% 
higher than that of SOLIS/VSM. This is due to the fact that 
HMI data includes small scale magnetic fields measurements. 
We study the magnetic pressure, p,„, in a longitudinal cross- 
section at about Q - 20° as shown in Figure [4] The overall 
pattern of appears to be the same when calculated from the 
HMI and VSM NLFFF model volume. In this study we found 
that the magnetic pressure of NLFFF model field from HMI is 
greater than that of VSM for same locations in the cross section. 
This is expected, since the magnetic pressure is proportional to 
the magnetic energy density as the magnetic energy of NLFFF 
model field from HMI is larger than that of VSM. The surface 
radial (Jr) and transverse (/,) electric current densities of the 
NLFFF field models based on HMI and VSM data are shown 
in in Figure |5] The value of the total radial surface electric cur- 
rent density flux of the NLFFF field models based on HMI is 
greater that than that of VSM. It agrees with fact that the HMI 
instrument measures more transverse magnetic field than that of 



VSM instrument. The transverse surface electric current density 
of the NLFFF field model based on HMI spreads more around 
the active regions than that of VSM as shown in Figure |5] This 
could reflect the fact ( |Pietarila et al.|[2012[ see) that the scal- 
ing factor between SOLIS/VSM and SDO/HMI is different for 
weak and strong fluxes. This difference is scaling factor may act 
as a weighting function when comparing electric currents de- 
rived from two models. In addition, the vector correlations of 
the radial and transverse surface electric current densities of the 
NLFFF field models based on HMI and VSM are 0.96 and 0.88, 
respectively. This indicate that there is more pronounced dis- 
crepancy in transverse electric current densities than radial one. 

5. Conclusion and outlook 

We have investigated the coronal magnetic field associated with 
full solar disk on November 09 2011 by analysing SDO/HMI 
and SOLIS/VSM data. We carried out nonlinear force-free coro- 
nal field extrapolations of a full disk magnetograms. The vec- 
tor magnetogram is almost perfectly flux balanced and the field 
of view was large enough to cover all the weak field surround- 
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Fig. 4. Magnetic pressure p,,, in the longitudinal cross-section at 
e = 20° for the (a) VSM and (b) HMI. The vertical and horizon- 
tal axes show radial distance in solar radius and longitude, (p on 
the photosphere, respectively. 



ing the active regions. Both conditions are necessary in order to 
carry out meaningful force-free computations. We have used the 
optimization method for the reconstruction of nonlinear force- 
free coronal magnetic fields in spherical geometry (Wiegelmann 
2007 Tadesse et al.|2009) l to campare the final NLFFF model 
field solution from HMI and VSM full disk data. 

We have found that the optimal value for smoothing param- 
eter is ii\ - 0.05 for full-disk HMI data for the purpose of pre- 
processing. We conclude that the choice v = 0.001 is the optimal 
choices for HMI full disk data set for our new code as investi- 

2012) 1. 



gated in [Wiegelmann et ak 



The magnetic field lines obtained from nonlinear force-free 
extrapolation based on HMI and VSM data have good corre- 
lations. However, the models show some disagreement on the 
estimated relative free magnetic energy which can be released 
during explosive events and surface electric current density. 
Reconstructed magnetic field based on SDO/HMI data have 
more contents of total magnetic energy, free magnetic energy, 
the longitudinal distribution of the magnetic pressure and surface 
electric current density compared to SOLIS/VSM data. Since 
the disagreement in free energy can be attributed to presence 
of weaker transverse fields in SDO/HMI measurements, it is not 
clear how important is the found (14.4%) diff'erence in free mag- 



140 



160 lao 

(a) SOL[S.'VSM 




ISO 



Fig. 5. Vector plot of the radial component of electric current 
density and vector field plot of transverse component of electric 
current density with white arrows. The colour coding shows 7, 
on the photosphere. The vertical and horizontal axes show lati- 
tude, Q and longitude, <p on the photosphere respectively. 



netic energy for flare and CME processes originating in magnetic 
fields higher in the corona. This aspect deserves a separate study. 
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